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ABSTRACT 

The relation between the zero-age main sequence mass of a star and its white- 
dwarf remnant (the initial-final mass relation) is a powerful tool for exploration of 
mass loss processes during stellar evolution. We present an empirical derivation 
of the initial-final mass relation based on spectroscopic analysis of seven massive 
white dwarfs in NGC 2168 (M35). Using an internally consistent data set, we 
show that the resultant white dwarf mass increases monotonically with progenitor 
mass for masses greater than 4M , one of the first open clusters to show this 
trend. We also find two massive white dwarfs foreground to the cluster that 
are otherwise consistent with cluster membership. These white dwarfs can be 
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explained as former cluster members moving steadily away from the cluster at 
speeds of < 0.5 km/s since their formation and may provide the first direct 
evidence of the loss of white dwarfs from open clusters. Based on these data 
alone, we constrain the upper mass limit of WD progenitors to be > 5.8M at 
the 90% confidence level for a cluster age of 150 Myr. 

Subject headings: white dwarfs — open clusters and associations: individual 
(NGC 2168) 

1. Introduction 

White dwarfs (WDs) are the final state of stellar evolution for the vast majority of 
intermediate- and low- mass stars. The upper mass limit of WD progenitor stars, M crit , is also 
the lower mass limit of core-collapse supernova progenitors. The current best observational 
estimate for M crit comes from spectroscopic analysis of WDs in the open cluster NGC 2516 
(Koester & Reimers 1996) and is M crit 8M ±2M . Due to the steepness of the initial-mass 
function (IMF), however, this range results in a factor of p=j 2 uncertainty in the number of 
supernovae and duration of supernova-driven winds resulting from bursts of star formation. 
This uncertainty in turn has a large impact on understanding the star-formation rate in 
galaxies (e.g. Somerville & Primack 1999), the evolution of starbursts (e.g. Leitherer et al. 
1999), and the fate of low- mass dwarf galaxies at early times (e.g. Dekel & Silk 1986). 

The best observational constraints on M crit are obtained from studies of WD populations 
in open clusters with ages < 150 Myr (Williams 2002). NGC 2168 (M35) is one of the richest, 
compact and nearby open clusters in this age range, with age determinations ranging from 
~ 100 Myr (von Hippel et al. 2000) to ~ 200 Myr (Sung & Bessell 1999). The WD cooling 
sequence of NGC 2168 has been discussed often in the literature, with recent photometric 
analyses by von Hippel et al. (2002) and Kalirai et al. (2003). 

With spectroscopy of cluster WD candidates, it is possible to determine unambiguously 
if the objects are WDs and, for the bona fide WDs, to determine T e fj and log g, and to derive 
the cooling age (r coo ;), mass and luminosity. For those WDs with cooling ages smaller than 
the cluster age and distance modulus consistent with cluster membership, subtraction of the 
WD cooling age from the age of the open cluster results in the lifetime of the progenitor 
star, and stellar evolutionary models can then be used to determine the progenitor star mass. 
Reimers & Koester (1988) applied this algorithm to two NGC 2168 WDs candidates identified 
on photographic plates by Romanishin & Angel (1980); however, their signal-to-noise was 
too low to determine the WD properties precisely. 
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As part of our ongoing program to identify and spectroscopically analyze the WD pop- 
ulations of open clusters, the Lick-Arizona White Dwarf Survey (LAWDS, Williams et al. 
2004), we have obtained high signal-to-noise spectra of eight candidate massive WDs in NGC 
2168. Six or seven of the observed WDs are hot, high-mass WDs likely to be cluster members. 
This sample doubles the number of known open cluster WDs with high-mass progenitors. A 
detailed photometric and spectroscopic analysis of the entire WD population of this cluster 
will be presented in a later paper. In this Letter we present the constraints on M crit and the 
upper end of the initial-final mass relation for WDs based on the high-mass WDs already 
analyzed. 



2. Observations & Analysis 

UBV imaging of NGC 2168 was obtained in September 2002 with the Lick Observatory 
Shane 3-m reflector and the PFCam prime-focus imager; additional UBV imaging of a larger 
field center on the cluster was obtained in January 2004 with the KPNO 4m MOSAIC 
camera. PSF-fitting photometry was obtained using the DAOPHOT II program (Stetson 
1987). Candidate WDs were selected by their blue excess in UBV color space. Figure 1 
shows the color-magnitude diagram of all objects detected in U, B, and V across the entire 
MOSAIC field. Several very blue, faint objects are observed in the diagram; these are our 
candidate WDs. All four WD candidates of Romanishin & Angel (1980) are recovered, as 
is the WD candidate in von Hippel et al. (2002). Astrometry and photometry for the WD 
candidates in this Letter are presented in Table 1. 

Spectroscopic observations of selected WD candidates were obtained with the blue cam- 
era of the LRIS spectrograph on the Keck I 10-m telescope (Oke et al. 1995). A l"-wide 
longslit at parallactic angle was used with the 400 1 mm" _1 , 3400A blaze grism for a resulting 
spectroscopic resolution of ~ 6A. The spectra were extracted and a relative spectrophoto- 
metric calibration applied using standard IRAF 1 routines. 

The T c fj and log g were determined for each WD using simultaneous Balmer-line fitting 
(Bergeron, Saffer, & Liebert 1992). The model spectra are updated versions of those in 
Finley, Koester, & Basri (1997). The WD evolutionary models of Wood (1995) were used 
to calculate the mass (M WD ) and cooling age (r WD ) of each WD. A distance modulus to 
each WD was measured by comparing the observed ^-magnitude to the absolute magnitude 



1 IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the 
Association of Universities for Research in Astronomy, Inc., under cooperative agreement with the National 
Science Foundation 
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M v calculated from the best-fitting model atmosphere and the appropriate WD cooling 
model. Errors in the fits were determined empirically by adding the noise measured for 
each spectrum to the best-fitting model spectrum convolved with the instrumental response. 
These simulated spectra were fit by the same method; nine iterations were used to calculate 
the scatter in T e g, log g, Mwd and twd- The fitting procedure is discussed in depth in an 
upcoming paper on the open clusters NGC 6633 and NGC 7063 (Williams et al. 2004). The 
atmospheric fits and derived WD masses and ages are given in Table 2. The Balmer line fits 
are shown in Figure 2. 

A systematic error in the fits of the hottest (T e ff > 50000K) WDs became apparent, as 
it was not possible to simultaneously fit all the Balmer lines. In these cases, the fits were 
limited to H/3, H7, and H<5. For LAWDS 22, no satisfactory convergence was achieved with 
these limited fits; the best-fit models are used in the analysis. WDs with T eS > 50000K are 
known to exhibit metals in the atmosphere and NLTE effects (e.g. Napiwotzki 1992; Holberg, 
Barstow, & Sion 1998), neither of which are included in the models. 

The progenitor mass for each WD was calculated by subtracting the WD cooling age 
from the cluster age. The age difference is the total lifetime of the progenitor star. The life- 
times of stars as a function of mass and metallicity are calculated from the stellar isochrones 
of Girardi et al. (2002). The progenitor mass for each WD likely to be a cluster member 
is given in Table 3 for an assumed cluster age of 150 Myr and for three different stellar 
evolutionary models: Z = 0.008 and Z = 0.019, both with modest convective overshoot, 
and Z = 0.019 without convective overshoot. Upper and lower errors are for la differences 
in twd- 

3. Discussion 

All eight of the observed WDs are much more massive than the typical WD mass of 
O.56M (Bergeron, Saffer, & Liebert 1992), and five of the WDs have apparent distance 
moduli (m — M)y « 10.5. It is therefore reasonable to assume that at least five of these 
objects are members of NGC 2168. LAWDS 11 is almost certainly not a cluster member. 
Its age is likely older than that of the cluster as a whole, and while the uncertainties in 
the spectral fits leave open the possibility that it is younger (which would require it to be 
hotter), the observed colors are more consistent with the cooler (and older) interpretation. 
The distance modulus of this WD is foreground to the cluster by little more than la, but the 
cooler temperature favored by the B — V color again favors the foreground interpretation. 

LAWDS 15 has a mass and age consistent with cluster membership, but the calculated 
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distance modulus is inconsistent with that of the other WDs by ~ 4a Mv - Assuming that the 
distance modulus is correct and that there is no difference in reddening between LAWDS 
f5 and NGC 2168, LAWDS 15 is ~ 185 pc closer than the cluster. Based on the spectral 
fits, LAWDS 15 has a cooling age of ~ 50 Myr. If the WD has been moving away from the 
cluster since its formation at a steady rate of only 0.4 km/s, it will have covered this distance. 
Therefore, it is possible that LAWDS 15 was once a cluster member and has escaped the 
cluster. 

The likelihood that a massive, hot WD would be found foreground along the line of 
sight to the cluster can be estimated from the luminosity function in Figure 16 of Liebert, 
Bergeron, & Holberg (2004). The luminosity function gives the space density of WDs with 
M > O.8M and r cool < 100 Myr as - 10~ 5 - 3 pc~ 3 0.5 mag" 1 . This results in an estimated 
0.1 hot, massive WDs in the ~ 30' x 30' MOSAIC field to a distance of 1 kpc. Therefore it is 
unlikely, but not impossible, that LAWDS 15 is a field WD. Based on these arguments, we 
will consider LAWDS 15 to be a cluster WD for this discussion. For similar arguments, we 
retain LAWDS 6 (~ 2a Mv closer than the other cluster WDs) as a likely cluster member. 

Figure 3 shows the initial-final mass relation of these seven cluster members, along with 
that of WDs from Hyades, Praesepe and the Pleiades (Claver et al. 2001) and from NGC 
2516 (Koester & Reimers 1996). Also shown are theoretical and semi-empirical data from 
plots in Claver et al. (2001) and sources therein. For the sake of consistency, the initial and 
final masses of each WD from the literature have been re-determined using our WD models 
and the published T eS and log g. 

From the figure, it can be seen that the NGC 2168 WDs form a monotonic sequence of 
more massive WDs originating from more massive progenitors. This conclusion is robust, as 
changes in the assumed age of NGC 2168 do not affect the relative positions of the points, 
only the absolute initial masses. Figure 3 also shows that the NGC 2168 initial-final mass 
relation agrees with the empirical relation derived from previous clusters. This need not be 
expected a priori, as NGC 2168 has a significantly lower metallicity than the other clusters 
in the diagram. Model core masses decrease with increasing metallicity, and the efficiency 
of mass loss processes could change with metallicity. 

LAWDS 15 and LAWDS 6, the potential escaped cluster WDs discussed above, fit 
the observed initial-final mass relation, providing additional evidence that they are cluster 
members. If these WDs are indeed escaped cluster members, they are crucial pieces of 
evidence that WDs can receive velocity kicks during mass loss and perhaps explain the 
observed deficit of WDs in other open clusters (see Williams 2004, and references therein). 
This and other potential explanations for these objects (e.g. binarity) will be discussed more 
fully in the later paper on WDs in NGC 2168. 
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The WD in Fig. 3 with the apparently low final mass is LAWDS 22. As mentioned 
above and visible in Fig. 2, the spectral fitting did not converge satisfactorily, despite the 
high signal-to-noise of the spectrum. This star is a close visual double with a redder com- 
panion (V = 19.12, B-V = 1.27, U-B = 0.93) 2"to the north. While resolved, this double 
is close enough that the spectrum of LAWDS 22 is likely contaminated by light from the 
neighbor star, resulting in the unsatisfactory fit. Light from the neighboring star may also 
be contaminating the photometric colors of LAWDS 22, which would also explain why the 
star lies redward of the 1M cooling track in Fig. 1. Given its measured T eS , LAWDS 22 
likely suffers from an extension of the systematic issue in the high-T e ff WD spectral fits de- 
scribed above. Other explanations for the location of this point could include magnetic fields 
(although no splitting is observed) or a low-mass, unresolved companion, but contamination 
by the neighboring star seems to be the most likely cause. 

Based on the NGC 2168 data alone, it is possible to place lower limits on the value 
of M cr it- Making the assumption that errors in the WD ages in Table 2 are Gaussian, we 
calculate that the oldest cluster WD ages are logr coo i > 7.76 with 90% confidence. For 
logTd = 8.15 and Z = 0.008, this corresponds to M crit > 5.81M at a 90% confidence level. 
This value is in agreement with that obtained by Koester & Reimers (1996). 

Improved constraints on the initial-final mass relation, including its intrinsic scatter and 
metallicity dependence, require improvements on existing observations. First and foremost, 
the ages of open clusters such as NGC 2168 must be determined to higher precision. Alterna- 
tives to main-sequence fitting such as lithium depletion studies or activity /rot at ion studies 
may provide the necessary constraints. Second, large samples of WDs from individual open 
clusters are needed to reduce the effect of systematics (such as errors in assumed ages) that 
plague the comparison of inter-cluster samples. The WD sample presented here is a start 
to that end, and the initial-final mass relation derived from these stars alone provides dra- 
matic confirmation of the existence of an initial-final mass relation, an idea that was strongly 
hinted at from previous open cluster studies and from other theoretical and observational 
work. Assuming that the majority of the remaining massive NGC 2168 WD candidates are 
cluster WDs, planned spectroscopic observations of these objects will soon result in a sample 
of nearly a dozen of WDs originating from a single stellar population, permitting, for the 
first time, studies of the intrinsic dispersion of the initial-final mass relation at high masses. 

M.B. and K.W. are grateful for support for this project in the form of the National Sci- 
ence Foundation AST-0307492. Any opinions, findings, and conclusions or recommendations 
expressed in this material are those of the author(s) and do not necessarily reflect the views 
of the National Science Foundation. 
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Table 1. WDs in the field of NGC 2168, photometry. 



LAWDS 


McCook & Sion 


RA 


Dec 


V 


ys y 


B-V 


ts—v 


Previous 


ID 


Designation 3, 


(J2000) 


(J2000) 








References 


NGC 2168:LAWDS 1 


WD J0608+242 


6:08:38.79 


24:15:06.9 


20.989 


0.019 


-0.035 


0.028 


1 


NGC 2168:LAWDS 2 


WD J0608+241 


6:08:42.30 


24:10:17.7 


21.569 


0.032 


0.061 


0.044 




NGC 2168:LAWDS 5 


WD J0609+244.1 


6:09:11.54 


24:27:20.9 


20.065 


0.017 


-0.128 


0.024 


1,2 


NGC 2168:LAWDS 6 


WD J0609+244.2 


6:09:23.48 


24:27:22.0 


19.863 


0.016 


-0.128 


0.023 


1,2 


NGC 2168:LAWDS 11 


WD J0609+241 


6:09:42.79 


24:11:05.4 


21.198 


0.025 


0.110 


0.037 




NGC 2168:LAWDS 15 


WD J0609+240 


6:09:11.63 


24:02:38.5 


20.785 


0.022 


-0.039 


0.032 




NGC 2168:LAWDS 22 


WD J0608+245 


6:08:24.65 


24:33:47.6 


19.657 


0.016 


0.008 


0.023 




NGC 2168:LAWDS 27 


WD J0609+243 


6:09:06.26 


24:19:25.3 


21.398 


0.026 


0.090 


0.039 


3 



a Format as in McCook & Sion (1999) I 

I 

References. — (1) Romanishin & Angel (1980), (2) Reimers & Koester (1988), (3) von Hippel et al. (2002) 



Table 2. WDs in the field of NGC 2168, spectral fits. 



ID 




dT cS 


log g 


dlogg 


M WD 


dM WD 


log(rwD) 


dlog(rwD) 


M v 


dM v 


(m - M) v 


LAWDS 1 


32400 


512 


8.40 


0.125 


0.888 


0.075 


7.386 


0.208 


10.454 


0.215 


10.535 


LAWDS 2 


31700 


1800 


8.74 


0.191 


1.072 


0.102 


7.929 


0.291 


11.132 


0.328 


10.437 


LAWDS 5 


52600 


1160 


8.24 


0.095 


0.824 


0.051 


6.158 


0.025 


9.542 


0.184 


10.523 


LAWDS 6 


55200 


897 


8.28 


0.065 


0.851 


0.036 


6.094 


0.027 


9.564 


0.108 


10.299 


LAWDS 11 


19900 


2792 


8.48 


0.367 


0.921 


0.213 


8.288 


0.463 


11.534 


0.769 


9.664 


LAWDS 15 


29900 


318 


8.48 


0.060 


0.934 


0.037 


7.693 


0.089 


10.754 


0.111 


10.031 


LAWDS 22 


54400 


1203 


8.04 


0.121 


0.721 


0.060 


6.169 


0.054 


9.155 


0.219 


10.502 


LAWDS 27 


30900 


500 


8.58 


0.164 


0.995 


0.086 


7.760 


0.220 


10.866 


0.310 


10.532 
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Table 3. Progenitor Masses for NGC 2168 WDs. 

ID Z = 0.008 Z = 0.019 Z = 0.019 a 

(Me) (M Q ) (M Q ) 

LAWDS 1 5.001°,;^ 5.08lg;?§ 4.86±g;?| 

LAWDS 2 7.02l?° 60 6.96l?° 49 6.63l?° 40 

LAWDS 5 4.628±g-ggi 4.734±g-°g? 4.540l°-°°l 

LAWDS 6 4.625+°;[™ 4.732l°;°°J 4.537l°;°° - 

LAWDS 15 5.5712m 5.62l°^ 5.371°,!° 

LAWDS 22 4.629l°;°°J 4.735l°;°°i 4.540l°;°°i 

LAWDS 27 5.8llJ;j£ 5.84lJ;j£ 5.591^ 



a No convective overshoot 







0.5 1 
B-V 



1.5 



2 



Fig. 1. — B — V, V color-magnitude diagram for NGC 2168. Squares indicate WDs presented 
in this study, both cluster members (filled) and the non-member (open). Solid curves indicate 
cooling curves for WDs with masses of O.4M (top curve), O.7M , and l.OM (bottom curve) 
at the distance and reddening of NGC 2168. The dotted line indicates the location of WDs 
with logr coo ; = 8.15, the assumed age of NGC 2168. 
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Fig. 2. — Balmer line fits for WD candidates in NGC 2168. The figures show the Balmer 
lines from H/3 (bottom) to H9 (top), with the curves showing the best-fit models. 
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Fig. 3. — Initial-final mass relation. Points with error bars are WDs from this study; other 
points are from the Hyades and Praesepe (Claver et al. 2001, open triangles), NGC 2516 
(Koester 1996, open squares) and the Pleiad WD (Claver et al. 2001, upper limit). The 
starred point with error bars is LAWDS 22, discussed in the text. Curves represent the 
theoretical initial-final mass relation from Girardi et al. (2000) (solid), the core mass at the 
first thermal pulse from the same models (dashed), and the quasi-empirical relation from 
Weidemann (2000) (dotted line). Slight horizontal offsets have been applied to the three 
WDs with Mi n it ~ 4.6M for the sake of clarity. Statistical error bars in the initial masses 
for these three WDs are smaller than the points. 



